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Abstract. We study the Klein paradox for the semi-classical Dirac operator on R with 
potentials having constant limits, not necessarily the same at infinity. Using the complex 
WKB method, the time-independent scattering theory in terms of incoming and outgoing 
plane wave solutions is established. The corresponding scattering matrix is unitary. We 
obtain an asymptotic expansion, with respect to the semi-classical parameter h, of the scat- 
tering matrix in the cases of the Klein paradox, the total transmission and the total reflection. 
Finally, we treat the scattering problem in the zero mass case. 
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1. Introduction 

In mathematics and physics, the scattering theory is a framework for studying and under- 
standing the scattering of waves and particles. The scattering matrix for the one-dimensional 
Dirac operator H is closely related to the transition probability of particles through a poten- 
tial. However, if the potential does not vanish at infinity, a Klein paradox might occur. The 
latter is of great historical importance in order to justify the existence of the antiparticle of an 
electron (the positron) and explaining qualitatively the pair creation process in the collision 
of particle beam with strongly repulsive electric field. The explanation of this Klein paradox 
usually resorted to the concept of "hole" in the "negative-energy electron sea". For more phys- 
ical interpretations we refer to Klein [14], Sauter (2Tj, Bjorken-Drell [2], Sakurai [20], Thaller 
[23] and Calogeracos-Dombey [3] for the history of the Klein paradox. This paradox appears 
also for the Klein-Gordon equation, here no concept of "hole" is needed (see Winter ^5] and 
Ni-Zhou-Yan [16j for a constant potential at infinity and Bachelot pLj for an electrostatic po- 
tential having different asymptotics at infinity). The comparison between the Klein paradox 
for this two equations has been discussed in [25l Part C]. A Klein paradox phenomenon occurs 
also in quantum field theory (see Hund [13] and Manogue [E]). It is clear that this paradox 
cannot appear for Schrodinger operators. 

For a scalar potential having real limits at ±00, the Klein paradox of the Dirac equation 
occurs if — > 2mc^. In this case the higher part of cr{H) intersects its lower part. If the 
energy E is in this intersection, for a wave-packet which comes from the left and moves towards 
the potential, a part of it is reflected, another part being transmitted. The transmitted part 
moves to the right and behaves like a solution with negative energy. Ruijsenaars-Bongaarts 
[To] (see also Thaller [23]) have mathematically treated the Klein paradox and the scattering 
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theory for the Dirac equation with one-dimensional potentials constant outside a compact 
set. They have established the connection between time-dependent and time-independent 
scattering theory in terms of incoming and outgoing plane wave solutions. The exact calculus 
of the scattering matrix for one-dimensional Dirac operator is only known for a few number 
of explicit potentials (see Klein [H] for a rectangular step potential and Fliigge [6] for the 
potential V{x) = tanh(2;)). Nevertheless, we are neither aware of works dealing with the 
asymptotic expansion of the scattering matrix, with respect to the semi-classical parameter h. 

For one-dimensional Schrodinger operators, there are several approaches which have been 
developed dealing with the computation of the transmission coefficient through a barrier. 
Ecalle [5j and Voros [24] have developed the so-called complex WKB analysis which gives 
approximations in the complex plane of the solutions of a Schrodinger equation. This approach 
is used in a new formalism by Grigis for the Hill's equation [llj. This method is also used 
by C. Gerard-Grigis [10] to calculate the eigenvalues near a potential barrier and by Ramond 
[T8] for scattering problems. For references and a historical discussion, we refer to Ramond 
|18| . The complex WKB method has been extended to a class of Schrodinger systems by 
Fujiie-Lasser-Nedelec [9]. 

The purpose of this paper is to give an asymptotic expansion, with respect to the semi- 
classical parameter h, of the coefHcients of the scattering matrix for the one-dimensional 
Dirac operator with potentials having different limits at infinity. We establish the exponential 
decay of the transmission coefficient in the Klein paradox case (cf. Theorem 12.11 below). We 
calculate the coefficients of the scattering matrix in terms of incoming and outgoing solutions. 
Therefore, we use the complex WKB analysis to construct solutions of the Dirac equation. 
The usefulness of this analysis is that it provides, rather than approximate solutions with error 
bounds, solutions in the complex plane with a complete asymptotic expansion with respect to 
the semi-classical parameter h, with a priory estimates on the coefficients. 

The paper is organized as follows. In the next section, we introduce the perturbed Dirac 
operator on M, study the time-independent scattering theory and state our main results. In 
Section [U we develop the complex WKB method and show a complete asymptotic expansion 
of the coefficients. In Section HJ the existence of incoming and outgoing Jost solutions is 
proved. In Section [5j we analyze the semi-classical behavior of the scattering matrix in the 
Klein paradox case. The total transmission over a potential barrier and the total reflection 
are studied in Section [6] and Section [71 Finally, in Section [HI we study the Klein paradox for 
the zero mass case. 

2. Assumptions and results 

We consider the self-adjoint Hamiltonian H = Hq + V, where Hq is the semi-classical Dirac 
operator on M: 

(2.1) Hq = —ihca h mc^P, 

ax 

with domain D{Ho) = H^{R) C Ti = L'^{R) C^, where /i \ is the semi-classical 

parameter, m > is the mass of the Dirac particle and c is the celerity of the light. The 
coefficients a, (3 are the 2x2 Pauli matrices satisfying the anti-commutation relation 

a/3 -|- Pa = 0, 

and = 0^ = I2, where I2 is the 2x2 identity matrix. 
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The operator V is the multiphcation by VI2, where F is a smooth electrostatic potential 
satisfying: 

(A): The function V is real on the real axis, analytic in the sector 

S = {x G C, |Im x\ < e|Re x| + ??}, 
for some €,rj > 0, and satisfies the following estimates: 
(2.2) \{V{x) - V^)\ = 0{{xy^) for Re (x) — > ±00 in S. 

Here, {x) = (1 + \x\'^)K ^ > 1 and V~ < . 



V- 



Fig. 1. The potential V 

The spectrum of the free Dirac operator Hq is ] — 00, —mc^] U [mc^, +00 [ and it is purely 
absolutely continuous. Under assumption (A) the operator H = Hq + V is sl self-adjoint 
operator and has essential spectrum (see Appendix lAl): 

(2.3) (TessiH) = ]- 00, -mc^ + y+] U [mc^ + y",+oo[. 

There are several representations of the matrices a, (3. For example, Hiller [12] used a = 
(72, 13 = (T3, Nogami and Toyoma [T^ used a = ct2, j3 = ai, where aj, j = 1,2,3, are the 
standard representation for Dirac-Pauli matrices. Most calculations with Dirac matrices can 
be done without referring to a particular representation (see Thaller [23l Appendix lA]). Here, 
we choose the 1 + 1 dimensional representation of the Dirac matrices 

(2.4) a = cTi=(^° J^, /3 = a3=(^J 
The solutions of 

(2.5) Hu= (^-ihcai^ + mc^a3 + V{x)l2^u = Eu, E eR, 

should behave as x — > ±00 like 
a^iE, h) exp(+-3-(m2c^ - (F± - Eff^x) + a±(S, h) exp(--^(m2c^ - {V^ - Eff'^x). 

Here, the square root (•)2 is to be defined more precisely according to the sign of m^c^ - 
{V^-Ef. 

In the following, we will use these intervals on the E-axis: 
I. y+ + mc^ < E, 

II. max(mc^ + V~ , — m(?)< E < + mc^, 
HI. V- +mc'^ < E <V+ - mc? if y+ - > 2mc^ 
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IV. V~ — mc^ < E < min(y~ + mc^, — mc^)^ 
V. E <V- - mc^. 
If V+ -V- > 1ra<?, the different regions are represented in the following figure: 



E k 



Region I 




Region IV 
Region V 

Fig. 2. Different regions on the E-axis 
We study the semi-classical behavior of the scattering matrix for the different values of the 
energy E. Let us describe now the time-independent scattering problem briefly. For -E G I, 
III or V, the four Jost solutions cj?j^, lj^jj (see Theorem 14. ip are the solutions of (12. 5p which 
behave exactly as 



(2.6) 



(2.7) 



out 



exp{=F7-«'(£^ 
he 

exp{±-^^>(£; 
he 



A(E - y±) 

^A{E-V^)-^ 

A(E - y±) 
±A{E-V^)-^ 



as X 



as X 



iboo, 



iboo, 



4/ E+m(? 



E—mc 



and sgn(S) = |^ for ^ [ 



2 2i 

-me , me \ 



with <^{E) = sgn{E)VE^ -m^c^, A{E) 
Analogous definitions of Jost solutions can be found in the works of Ruijsenaars-Bongaarts 
[To] and Thaller [23] for one-dimensional step potentials. In this paper, we denote ^/x, ^fx 
the positive determination of x G — > (a^)^, (x)4 respectively. 

The ordinary scattering problem is the following: what are the components of a solution u 
of the Dirac equation (|2.5p in the basis (w^^ti "^out) of the outgoing Jost solutions, knowing 
its component in the basis (wj", w^^) of the incoming Jost solutions. The 2x2 matrix relating 
these coefficients is called the scattering matrix and we will denote it by 

S21 S22 



Precisely, if we take u a solution of (|2.5p . 
the scattering matrix is such that 



out ' 



t^out 
^out 



which is equivalent to 
(2.8) 



{^in^^tn) = (^c]ut'^out)S- 
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Since V is real on the real axis, we have (see (|2.4p ) 

(2.9) ^ = (3u;t^,. 

We also have the following relations between the coefficients of S{E, h): 

(2.10) sii{E,h)=S22{E,h) and si2{E,h) = -s^{E,hf:^^^^, 
so that sii and S12 determine completely the scattering matrix. 

The reflection and transmission coefficients R[E, h) and T{E, h) are, by definition, the 
square of the modulus of the coefficients S21 and su respectively. They correspond to the 
probability for a purely incoming-from-the-left particle to be reflected to the left or transmitted 
to the right. Using (12.91) . (I2.10p and calculating the determinant of (|2.8p . we have the well- 
known relation R{E, h) + T{E, h) = 1 and, the scattering matrix S{E, h) is unitary. 

To calculate the scattering matrix S(-E, h) we will use the transfer matrix T, which is defined 

by 

The determinant of this matrix is equal to 1 since the two Wronskians >V(a;j~, u;"^^) and 
>V(a;^^, (jj^ are equal to —2 (see Definition 13.41) . Using the relation (|2.9p . we obtain that T 
is determined by two coefficients: 

t{E,h) r{E,h) 
—{E, h) J{E, h) 

Moreover, using that det(T) = 1, we obtain 

(2.12) \t{E,h)\'^ -\r{E,h)\^ = 1. 

Consequently, we can write the scattering matrix in terms of the coefficients of the transfer 
matrix T: 

(2.13) S = =,^ ( ' 



(2.11) T 



T(S, h) \ -{E, h) 1 

We will use WKB approaches to describe the amplitude of the coefficients of the scattering 
matrix for h\0. For these, let us introduce the following definition. 

Definition 2.1. (See Sjdstrand {22\) A function f{z,h) defined in C/x]0,/io[, where U is an 
open set in C and ho > 0, is called a classical analytic symbol ( CAS) of order m G in h if 
f is an analytic function of z G U and if there exists a sequence {aj(z)) of analytic functions 
in U such that 

• For all compact set K C U, there exists C > such that, for all z & K, one has 

\a,{z)\<C^+'f. 

• The function f{z,h) admits the series X^o<j<+oo '^^ asymptotic expansion 
as h goes to zero in the following sense. For any Ci > C , we have 

f{z,h)- a^{zW^^ = 0{e-P'^), 

0<j<h-^/eCi 

for some p > and all z G K 
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The main theorem concerning the Klein paradox case for m > (i.e. for the energy level 
E G III) is the following: 




Fig. 3. Graph of V{x) + mc and V{x) — mc 



Theorem 2.1. [Klein paradox] Let V be a potential satisfying assumption (A), E G III and 
m > 0. Suppose that there are only two simple zeros ti{E) < t2{E) of m^c^ — {V{x) — E)'^ 
(see Fig. 3). Then there exists three classical analytic symbols 4'i{h), 4>2{h) and <p^{h) of 
non-negative order such that: 

(2.14) sii = S22 = {l + Hi{h))e^v{-S{E)/h}e^v{iT{E)/h}, 



(2.15) S2i = (i + /i02(/i))exp<! ^ Ui(E)VE- + 



rti{E) ^ 
/ Q-{t,E)dt 

J —OO 1 



(2.16) si2 = (^ + h4>^{h)) exp <! ^ I t2{E)VE+ 



r+oo 
/ 

Jt2{E) 



Q+{t,E)dt 



where S{E) is the classical action between the two turning points ti{E) and t2{E) 

S{E) 



*2(-B) /™2„4 _ 



ti{E) 



w?d^ - {V{t) - Ef 



dt. 



Moreover 



and 



where 



Q-{t,E) 
Q+{t,E) 

'ME) 



'-mV + (F(i) 


-Ef 




'-mV + (F(i) 


-EY 



for t < ti{E), 
for t > t2{E), 



/ti(i^) r+oo 
Q-{t,E)dt- / Q+{t,E)dt + ti{E)VE^ + t2{E)VE+, 



, -m^c^ + (y± - E)2 
E = . 
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We remark that this scattering matrix behaves like in the case of the Schrodinger operator 
with a barrier potential. In particular the term e~^^^^^'^ which decays exponentially, can be 
viewed as a tunneling effect (see Ramond [IBl Theorem 1]). 



In the zero mass case, we have the following theorem: 

Theorem 2.2. [Zero mass case] Let V be a potential satisfying assumption (A), E € III 
and m = 0. Suppose that there is only a simple zero tQ{E) of V{x) — E. Then, there is a 
classical analytic symbol 4>{h) such that: 

(2.17) sii = S22 = (l + /i'/'(/i))exp{iro(^)//i}, 



(2.18) S2i = Oih), 

(2.19) si2 = 0{h). 
Here, 

To{E) = T{E) 




-V{t) + V-)dt - / {V{t) - V+)dt + to{E){V+ -V-)]. 

Jto(E) J 

Remark 2.1. We can not permute the limits of the scattering matrix S as m ^ and h ^ 0. 
Indeed, if we take the limits of Si2 in ^2. 1 0^) and ^2.19^) . we obtain 

lim lim |si2| = 1, lim lim |si2| = 0. 

Now, we come back to the non-zero mass case and we treat reflection and transmission 
cases (see Sections [71 [6|) . 

If we take the energy level E GlI, there are two Jost solutions Wj", a;~^^ satisfying (12. 6|) and 
(12.71) for X — > — oo and there does not exist an oscillating solution for x — > +oo. Instead, 
as x ^ +00, there exists an exponentially decaying solution and an exponentially growing 
solution. Since the last function doesn't represent a physical state we limit ourself to the one 
dimensional space generated by the decaying solution wj" (unique up to a constant). This 
function satisfies (see Theorem |4JJ : 



(2.20).;+ ~ exp{-lVmV-(y+-i^)M J^^J^ I as x ^ +oo. 



41 mc^+E-V+ 



-E-V^ 



In this case we have 



Theorem 2.3. [Total reflection] Let V be a potential satisfying assumption (A), E ^11 and 
m > 0. Suppose that there is only a simple zero ti{E) of m'^c'^ — {V{x) — E)"^ . Then the vector 
space of the solutions of {H — E)u = with u bounded is a one dimensional space generated 
by 



U = ^u, + "out'^ouf with 
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(2.21) = -i(l + /i(/>i(/i))exp|^ ^ Q-(t, + V^ti(^)^ | . 

Moreover 

u = with 

( 1 / ftiiE) rh(E) ^ 

exp<^- / Qt{t,E)dt + V-E+ti{E)+i {t, E)dt + iVE^ti{E) 

Here 



Q^{t,E) and E^ are the functions of Theorem \2.1\ and (t)j{h), j = 1,2 are classical analytic 
symbols of non-negative order. 

For E S IV, there is also total reflection cases which can treated similarly to the previous 
theorem. As in |l9j. there is also a scattering interpretation of the previous theorem. Since 
we work in a one-dimensional space, the scattering matrix is now a scalar. 

Remark 2.2. [Scattering interpretation] We call = a;^ + a^ut^o^t "incoming" 
solution. In the same way, there exists a unique hounded solution 

which is called the "outgoing" solution. 

If u is a bounded solution of {H — E)u = (i.e. u = Au^^) then u = Bu°^^ . The scattering 
matrix S is defined by 

B = §A. 

Prom 112.21]} . we have 

S = «out = -^(1 + |x ( /'^""^ ^^"^^ + VE^t,{E)j I . 

For -E G I or V, a total transmission phenomena occur: 

Theorem 2.4. [Total transmission] Let V be a potential satisfying assumption (A), E £ 
I, m > and rn^c^ - {V{x) - Ef / 0. Then there are a classical analytic symbol (p{h) and 
positive constant C such that: 

(2.22) sii = S22 = {l + h(P{h))e^p{if{E)/h}, 

(2.23) S21 = 0{e-^l^) and si2 = 0{e-^'^), 
where 

Q-{t,E)dt+ / Q+{t,E)dt, 

-oo Jo 

and Q~{t,E), Q^{t,E) are the functions of Theorem \2.1\ defined here for any t E M. 

We can calculate the scattering matrix for -E S V in the same way of E £ I. We remark 
that the behavior of the incoming and outgoing Jost solutions exchanges between these two 
cases. This is in agreement with the physical interpretation (see [23l p. 121]). 
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3. Complex WKB solutions 

We wish to find a representation formula for the solutions of (|2.5p . from which it is possible 
to deduce the asymptotic expansion in h. The method is known as complex WKB method. 
See [18] [lO], [7], [8], [9] for constructions of solutions of the Schrodinger equation. 

In a complex domain 5, we study the Dirac system (|2.5p which is of the form 

-.X /TT / \ ( mc'^ + V(x) — E —ihc4- \ , ^ 

3.1 {H-E)u{x)={ 2 Trf\ ]u{x)=0, 

or equivalently 

<'-^) 7S"(^) - ( 'V ) "(^)- 

where v{x) = ( ^ q ) u{x) = M^^u, and the functions 



-mc^ =F (V(x) - E) 
9±{x) = 



c 

are holomorphic in S. The following considerations will lead to the construction of complex 
WKB solutions for Dirac system. 

3.1. Formal construction. First, we introduce a new complex coordinate 



(3.3) z{x) = z{x,xo)= {g+{t)g-{t))-2dt= {g+{t)g^{t))-2dt, xoeD. 

Jj{xo,x) Jxo 

One of our tasks will be of course to choose the simply connected subset D of 5 such that 
t — > {g-^-{t)g-{t))2 is well-defined, but let's work formally for a while. The 7(xo,x) is any 
path in D beginning at xq and ending at x. 



Definition 3.1. The zeros of the function 

■rr?c^ - {y{x) - Ef 



g+{x)g-ix) 



c2 



are called the turning points of the system II 3. 
Definition 3.2. For x fixed in D, the set 

S^yeD, Re j\g+{t)g^{t))^2dt = 

is called the Stokes line passing through x. 

_|__2_ 

We look for solutions of the form e hw±{z). We note that due to the possible presence of 
such turning points, the square root in the definition of z{x) might be defined only locally. 
By formal calculations, the amplitude vector w±{z) has to satisfy 

(3.4) ld-z'"^^'^ = [-Hizr ±^ 
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The function H{z{x)) is given by 

(3.5) ^— ^+(>'o.o-^)^'" 



g+{x) J \—mc^ — {V{x) — E) 

for z{x) in an open simply-connected domain of the z-plane, where z — > H{z) is well-defined 
and analytic. 

In order to obtain a decomposition with respect to image and kernel of the previous system, 
we conjugate by 

P (,^-\( H{z) ±iH{z)-^ \ 1 _ / H{z)-^ H{z)-^ 

^^^> ~ 2 \ H{z) TiH{z)-'^ i ' ^± \ ^iH{z) ±iH{z) 



and obtain a system for w±{z) = P±{z)'w±{z), 



d _ / 



H'{z) 



(3.6) i^±(^)=l H'i.) ^ ]^±i^l 
where H'{z) is shorthand for ■^H{z). The series ansatz 

(3.7) w^iz) = y^( ^^n,±{f \^ 



with Wo,± = 1 and, for n > 1, the recurrence equations 
/ X f d 2\ , , H'(z) 

(3.8) [ — ± - \ W2n+l,±{z) = rj, . W2n,±{z), 



dz hj ' H{z) 

^, d , ^ H'iz) 

(3-9) 3-W2n+2,±(^;) = . W2n+IMZ), 

dz "\z) 
give us a formal solution up to some additive constants. The solutions are fixed by setting 

Wn,±{z) =0, n > 1, 

at a base point z = z{x) where x G is not a turning point. We note that the previous 
equations for Wn,± are similar to the ones obtained by a complex WKB construction for scalar 
Schrodinger equations. See for example the works of C. Gerard and Grigis [10] or Ramond 

m 

Let O = J7(£') be a simply connected subset of D which does not contain any turning point. 
Then the function z = z{x) is conformal from VL onto z{Q.). Assume that z G z{Q). If T±{z, z) 
denotes a path of finite length in z{VL) connecting 'z and z G z{Vl), we can formally rewrite the 
above differential equations for n > as 

W2n+l,±{z) = / exp{±-{C - z))——W2n,±{C)dC, 

Jr±{z,z) ^iU 

W2n+2,±{z) = / -E777T'^2n+l,±{0dC, 

Jr±{z,z) -"lU 
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or after iterated integrations, as 



W2n+l,±{z) = / ■ ■ ■ / *5^P ( =tr(Cl - C2 H h C2„+l - Z)] X 

JT±(z,z)Jv±(z,C2n+l) -'r±(2,C2) V / 

^'(Cl) ^'(C2n+l) 



dCl • • • dC2n+l, 

I) 

W2n+2,±{z) = I I ■ I exp { ±^(Cl - C2 H C2n+2) ) X 



r±{2-,z) 7r±{2,C2n+2) -'r±{i-,C2) ^ ^ 



i7(Cl) i^(C2n+l) 

/r±{i-,C2) V 

g^(Cl) g^(C2n+2: 
H{Cl) ■ ■ ■ ^(C2n+2) 



dCl • • • dQ2n+2 ■ 



3.2. Convergence, /i-dependence and Wronskians. We now give to the preceding formal 
construction some mathematical meaning in simply connected, turning point-free compact sets 

Lemma 3.1. For any fixed h> 0, the series \3. 7]) converges uniformly in any compact subset 
of and 

(3.10) <^°(X, h)=Y, W2n,±{4x)), wf'^ix, h)=Y, W2n+l,±{z{x)), 

n>0 n>0 

are holomorphic functions in 1^. 

Proof. By assumption on Q and on V, the functions Wn,± are well-defined analytic functions 
in Q. For compact subsets K C Q and z , z G z{K) there exist positive constants C^{K) > 0, 
depending on the semi-classical parameter h and the compact K such that 

(3.11) sup exp(±^C)S <CiiK). 



If we denote the maximal length of the paths T±{z, ■) C K in the preceding iterated integra- 
tions by 

L = _ max inin \^(z, z)\ < 00, 

z,zG z{K) 7(2,2) 

then 

C^(i^)"L" 



sup \wn,±{z)\< j , n > 0, 

z&z{K) 



where the bound ^ comes from the volume of a simplex with length L. Then, the lemma 
follows. □ 

Thus, we have the uniform convergence of the series (|3.7p for w±{z) and complex solutions 
(3.12) u^{x) = e^^T^izix)) ^|dd(^) 

of the original problem (|3.1|) on any turning point-free set where 

(3,13, ^ p^w -^yv 
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We write these solutions u±{x) as 

(3.14) u±{x;xo,x), 

indicating the particular choice of the phase base point xq, in (|3.3I1 . which defines the phase 
function z{x) = z(x;xo), and the choice of the amplitude base point z = z{x), which is the 
initial point of the path T±{z, •). 

Definition 3.3. For x £ Q, fixed, we define 0± = Q±(5?) the set of all x G O such that there 
exists a path T±(z{x), z{x)) along which x — > ibRe z{x) increases strictly. 

Proposition 3.1. The functions Wn,± are classical analytic symbols of order [^^] in 0±. 
The functions w^^'"'{x,h) and Wj^'^{x,h) given by the identities L3.10\} are classical analytic 
symbols of order and 1 respectively in Moreover, we have for any N £ 'N, 

N 

w'^''''{x,h)-^W2n,±{z{x)) = Oih^'^'), 
n=0 

N 

W^^\x,h)-^W2n+l,±{^{x)) = 0(/i^+2), 
n=0 

uniformly in any compact subsets of fl± . In particular, 

w^'^'^ix, h) = l + Hih), w°^\x, h) = h^{h). 

Here and in all this paper, 4>{h) is a classical analytic symbol of non-negative order not neces- 
sarily the same in each expression. 

The proof is just the same as that of [TOl Prop. 1.2] and [H Prop. 3.3]. The key point is the 
following: since the iterated integrations defining Wn,±{z) contain terms of the form exp(ib^), 
one has to make sure that x i — > itRe {z{x)) is a strictly increasing function along the path 
r±(2^)'2). In other words, the paths T ±{z{x) z{x)) have to intersect the Stokes lines, that is 
the level curves of x i — > Re [z{x))., transversally in a suitable direction. 

Definition 3.4. One defines the Wronskian of two C'^-valued functions u = {ui,U2), v = 
{vi,V2) as 

W(n, v) = U1V2 — U2V1. 

Remark 3.1. For two solutions u and v of the equation h2. 5|) . the Wronskian W(ti,f) doesn't 
depend on x and is zero if and only if u and v are proportional. 

If w = au + (3v with a, [3 £ C, then 

_ W{w,v) _ W{w,u) 

Elementary computations give the following complex Wronskian formulas for complex WKB 
solutions with different phase and amplitude base points in terms of w^^"^ and 

Lemma 3.2. Let xq and yo be two points m = Q{E). If, for given x and y, the canonical 
sets Q±{x) and ^±{y) have a non-empty intersection, then for any x E ^l±{x) n ^±{y) one 



KLEIN PARADOX AND SCATTERING THEORY 



13 



has 

(3.15) W{u±{x;xq,x) ,u±{x;yo,y)) = ±2iexp (^±^{z{x;xo) + z{x;yo)) 

X (w5^^^(x;xo,x)u;|''™(x;yo,y) - w''^''''{x;xo,x)w'^'^{x;yo,y)) ■ 

If, for given x and y the canonical sets Q±{x) and ^jj^iy) have a non-empty intersection, then 
for any x £ il.±{x) n ^i=p(y) one has 

(3.16) W{u±{x;xo,x) ,u^{x;yo,y)) = ±2zexp (±^{z{x;xo) - z{x;yo)) 



h 

X (<<i<i(x;xo,^)<<^<^(x;2/o,y) - w^r°(^;^o,£)<'°(x;yo,y)) • 

4. JOST SOLUTIONS 

The Jost solutions of Hu = Eu, are characterized by the behavior of the solutions at 
infinity. We construct here the Jost solutions copying the procedure described in Section [31 
the new point here being that the solutions we seek are normalized at infinity. In all this 
section we will work in two unbounded, simply-connected domains ^1~{E), ^1~^{E), where 
Re {V{x) + mc^) < E, Re {V{x) — mc^) > E respectively and which coincide with S for 
Re X sufficiently large. The existence of such domains is of course an easy consequence of the 
behavior of V at infinity in S (see assumption (A)). 

First we define the phase functions with base point at infinity, 

(4.1) z(x,±oo)=/ - dt 



J itoo V / V / 

We also see that the integral converges absolutely, hence 

/^2^4_(y±_^)2x 1/2 

z{x, iboo) = ( 5 X + o(l), [x — > itoo). 



If the determination of the square root in z(-, •) are the same, we get the following equalities 

z{ti, ibcX)) =z{x, iboo) — z{x, tl) = z{t2, ±Oo) — z{t2, ti) 

(4.2) 

z(tl, — OO) — z{ti, +00) =z{x, —00) — z{x, +00) = z(t2, —00) — z{t2, +Oo), 

where z{-, •) is defined in (|3.3p . (14. 1|) and x, ti, t2 £ D. 

Next we define the amplitudes based at infinity. We will only define the amplitudes at +co 
since the situation is similar at —00. As in Section 3 of [18j, we choose infinite paths 7±(a;) 
starting from infinity and ending at rr, which are asymptotically like lines of the form {\mx = 
=F/oRex} for some /? > 0, such that x 1 — > ^'Kez{x) are strictly increasing functions along 
7±(x). Denoting the path z{'^±{x)) by T±{+oo^ z{x)) and setting wq^± = 1, we inductively 
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define Wn,±{z) by 



W2n+l,±iz) = / exp{±-{C - z))——W2n,±iOdC, 

Jr±{+oo,z) ^lU 

W2n+2,±iz) = / T777T^2n+l,±(C)c^C, n > 0. 

ir±(+oo,2) ^UJ 



Noticing that 



H'{x) _ mc^ V'{xj _ ^,,,^s 



H{x) 2 {V{x) - EY - m?c* 



0{{x) ), (5 > 1, as \x\ — > oo, 



one constructs well-defined complex WKB solutions u^^ corresponding to these base points, 
proceeding as in Section [3l Here, I and r stand for left and right and correspond respectively 
to X — > — oo and x +oo. Up to a constant pre- factor, u^i{x) are the previously defined Jost 
solutions: 

Lemma 4.1. Let u^i{x) be the complex WKB solutions with phase and amplitude base point 
at infinity. Then 

(4.3) utix) ~ exp(±-l(mV-(y+-i?)2)V2^)(^ ^^^^^ x ^ +oo, 

(4.4) nf(x) ~ exp(±^(mV-(y--i?)2)V2^)(^ ^ 



X — > — oo, 



with 



-y± + E - mc2 



Proof. We just check the asymptotic behavior of w^;(x) at infinity. Since H{z{x)) — > as 
X — > iboo, using (I3.12P and (|3.13l) we get by an elementary calculation 



t^{x) ~ exp(±-l(m=c' - (V+ - ( f ^l "+ ) ( J ) ■ - ^ +»• 



This ends the proof of lemma. □ 

Let us now choose the determinations of (m^c^ — {V^ — E)'^)^/'^ and according to the 
intervals on the i?-axis. This fixes the choice of uf^ and we can construct Wj^, satisfying 





±ia~^ 


l/a+ 


l/a+ 




zizia~ 


1/a- 


1/a- 




-Ef) 



M and (I2711 . 

L For G I, we choose (m^c^ - (y± - ^)2)V2 ^ G e*''/^M+ and we denote 

(4.5) u;,;:=e-/V, < := -e^/V, c^o+ut := ^^^^V, ^"ut := -e^^N"- 

2. For S G II, we choose (m2c^-(y--£;)2)i/2 g ^]^+^ ^j- ^ e^''/%+, {w?c^ -{V+ -Ef f/"^ G 
M~, a+ G M+ and we denote 

(4.6) <:=e-/V, <ut "^^^V, := n+. 

3. For E G III, we choose (m^c^ - (F± - £;)2)i/2 g G e=F*'^/^M+ and we denote 

(4.7) u;r:=e*^/V> -e^/V, u;o\t := e'^'^/V, u^'t := -e^'^V- 
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4. For E G IV, we choose {m^c^ - {V+ - EfY''^ G a+ G 6"*^/%+, (m^c^ - (F" - 
^)2)i/2 e R-, a+ e M+ and we denote 

(4.8) u;.+ :=-e— /V, ^0"^* := e'^'^/V, ^d--=nT- 

5. For G V, we choose {m^c^ - {V^ - E)^y/^ G iM+, a± G e-*''/''M+ and we denote 

(4.9) u:r := /^^^^ ^in a;+ , := e— /^n,", u;" , := -e— Z^^. 

Theorem 4.1. For rea/ E, h2. 5|) /tas solutions of the following form: 

1. For E £ I, III or V, there are four Jost solutions uj^, w^^- which behave like 

(4.10) u:i ~ exp{T-^$(i?-y±)x}f ^fjf "3 1 1 as x ^ ±00, 



(4.11) <ut ~ exp{±— $(^-y±)x} { ^4;p_,^±^i ) asx^ioo, 



A(F- 


y±) 


TA(i?- 


y±) 


A(F- 


y±) 


±A{E- 


y±) 



with ^E) = sgn{E)VE'^ -m'^c'^, A{E) = and sgn{E) = ^ for E ^ {-mc^ ,mc\ 

2. For G II (resp. -EG IV^, there are two Jost solutions lo^, uj^^^ (resp. lo^, w^^^-J which 
behave as in 14 -10^ , U-H^ and a decreasing solution uj'^ (resp. ^^)- The solutions oj^ behave 
exactly like 



(4.12) ~ exp{T-VmV-(y±-i?)M ^r^^ \ as x ^ ±00. 



vac- 



■^+E-V± 



For G I, III or V, according to the relation (|2.13|) . it is sufficient to calculate the two terms 
r{E, h), t{E, h) in T to obtain the matrix §. The definition of the Wronskian (see Definition 
[3T4II leads to: 



(4.13) tiE,h) 

(4.14) r{E,h) 



+ ) 



5. The Klein paradox case 

We suppose that V satisfies assumption (A), the energy E G III and m > (see Fig. 2). 
In this section we will work in two unbounded, simply-connected domains Q~{E), Q~^{E), 
where Re {V{x) + mc^) < E, Re {V{x) — mc^) > E respectively and which coincide with S 
for |Re x\ sufiiciently large. Using Theorem 14.11 Proposition 13.11 and (|4.7|) there are two Jost 
solutions in Q,^{E): 

(5.1) = exp{lz(x, ±00)} ( ) (1 + 

r("'^°")H±^(.(x)) 

The function z{x, ±00) is defined by (|4.ip and 



(5.2) u;^^, = exp{— z(x,±oo)}( , -^XH-i ]il + hm)- 



(5.3) H{zix)) 



E - V{x) + mc^X 4 



E — V{x) — mc 
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On n^{E)r\R, we have: 

(5.4) 4.,±») = .r M^-mr-'nviiE-v-r-m^c'^^ 

J ±00 V c V c 



(5.5) H{z{x)) 



I E -V{x) + mc^ 
E — V{x) — mc^ 



We suppose that there are only two real turning points ti{E) < t2{E) and that they are 
simple. Notice that ti{E) is a zero of E — V{t) — mc^ and t2{E) is a zero of E — V{t) + mc^. 
In that case the Stokes lines are as shown in Fig. 4. In order to obtain S, we compute the 
Wronskians given in (|4.13ll . (|4.14l) and then the coefficients t{E,h), r{E,h). 



71 


/ / 
















/t2{E) 


72^ 







Fig. 4. The turning points and the paths -yj 



Computation of W{u!'^^^^,uj^): Since the two solutions w^^jWj^ are defined in we 
can compute this Wronskian in Q~^(E) and from Lemma [3?2l we obtain 

(5.6) WK+t,^i+) = -2. 

Computation of W{uj^,lo^): The two solutions u>^,uj^ are defined in Q~ (E) , il.^ (E) re- 
spectively. Since the Wronkians W{uj^,uj^){x) are independent on x (see Remark l3.ip we 
compute this Wronskian in Q~{E) for example. For that we extend uj^ , which is defined in 
Q~^{E), into Q~{E). We will extend the square root in uj^ which is defined in ^~^{E), into 
C \ {-Im {z) > 0, Re {z) = ti{E)} U {-Im {z) > 0, Re (z) = t2{E)}. Thanks to the struc- 
ture of the Stokes lines between ti{E) and t2{E), we can find a path 71 from +00(1 + i6i) to 
—00(1 — i6i) (for 61 > 0) transverse to the Stokes lines along which we can extend lo^. We re- 
mark that between ti (E) and t2{E) on the real axis we have {E-V{t))^-m'^c^ < 0. The exten- 



Sion 0ftG]t2{E),+00[-^ -m c pQj^^i^g ^i^j^ -i m c -(^-VW) ]t^(^E) , t2iE)[ 



and with —\J^^—^^^^ — ■m_c_ j _ qq^i-^(^e)[. On the other hand, the extension of H{z{x)) 
stay in on ] — 00, ti (£')[. If we denote by ujf^^ the extension of uj^^ along 71, we have: 

^in'' = eM\{-z{xM{E)) + z{t2{E),+^)+S{E))}(^_^^^^^^^^^ 
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with z{t2{E),+oo) defined in (|5.4p and 



(5,7) z(x,ti(E)) = i/ I '-^ '"''jf ] dt. 



(5.8) s(E) = f'"' y^!ZZIpZM!,, 

Here, (^(^zXM)!^!!^) ^ g m+ for t g] - oo,ti(i?)[. 
Then, 

WK'^^i) = -2(1 + exp{i(z(ti(i^), -oo) + ^(t2(i?), +oo) + S{E))}, 



where z{t2{E),+oo), z{ti{E), —oo) are defined in (15.41) . 

Computation of W{u!~^^,io^): This wronskian is also between two solutions which are de- 
fined in different domains, then we extend one of these solutions into the domain of the other 
solution. For example we extend u!^, which is defined in Q^{E), into Q~{E) which is a subset 
of C \ {Im (z) > 0, Re (z) = ti{E)} U {-Im (z) > 0, Re (z) = t2(^)}. Here, we can also find 
a path 72 from +00(1 + 162) to — oo(l + ^(^2) for 62 > transverse to the Stokes lines along 
which we can extend into Q~{E). If we denote by uj^'^ the extension of u)^ along 72, we 
have: 

u^t/ = eMliz{xME)) + z{t2{E),+^) + Sm}l^.f^^^^^^^^^ ) (1 + 

Here H{z{x)) G M+ on ] -oo,ti(^)[ and z{x,ti{E)), S{E) are defined in ([52]), (iBlHIl . 
The computation of W(w~^^, W;^) yields: 



y^i^nt^^t) = 2^(1 + exp{-i-zit,iE), -00) + zME), +00) + SiE))}. 



1 

h 

Then, we obtain (see (l4T3]l and ([4141) 1 



tiE,h) = {l + h^{h))exp{^{z{ti{E),-oo) + z{t2{E),+oo)+SiE))}, 
r{E,h) = -i{l + h^ih))exp{^{-z{ti{E),-oo) + z{t2{E),+oo) + SiE))}. 
Since (p{h) is a classical analytic symbols of non-negative order and using (I2.13P we have: 
sii = =^^ = {l + hMh))exp{^iziti{E),-^)+zit2{E),+^)-S{E))}, 

S21 = ^^^ = {^ + hMh))exp{l{z{t^{E),-^)}, 

^ ~T{Eh) ^ (' + '^Mh))exp{^{z{t2{E),+^)}. 

The functions 4>i{h), (/>2(/i), 4>?,{h) are classical analytic symbols of non-negative order. This 
ends the proof of Theorem 12.11 
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6. Total transmission 

We suppose that V satisfies assumption (A), the energy E £ I oi E € V and m > (see 
Fig. 5, Fig. 6). 



E 




V + mc^ 
V~— mc^ 




Fig. 5. V ± mc^ and E Gl 



Fig. Q. V ± mc^ and E gV 



We suppose that there exists no real turning point. In that case the Stokes hues are hor- 
izontal lines near the real axis. We will only work for E £ I. The case where E £ \ can be 
treated similarly. 

In this section we work in Q~{E) defined in the previous section. Now this set is a neigh- 
borhood of the real axis. Using Theorem 14.11 Proposition 13.11 and (|4.5p there are four Jost 
solutions: 



(6.1) 



'^in =exp{— z(x,±CX))} 



-±1 



out 



exp{—z{x, ±oo)} 



h 



H{z{x)) 

H{z{x)) 
±H{z{x))-'^ 



{l + h4>{h)) 
(l + /i0(/i)). 



The functions z(x,±oo) and H{z{x)) are defined in (14.1(1 and (15. 3|) and coincide with (15. 4|) . 
(15. 5|) on the real axis. Here, the setting is different from the previous section. The solutions 
u!^ and are defined in the same domain i}~{E) and there are no problem to extend the 
different square roots. 

As in SectionlU it is sufficient to calculate the two terms r{E, h), t{E, h) (see (|4.13l) . (|4.14l) ) 
to obtain the matrix S. 

Computation o/yV(Wjjj, Wj^J, yV{u}'^^^^uJ^ : Since the function H in (|6.ip is the same for Wj" 
and Wj^, we have: 

(6.2) W«,a;+)(x) = -2(l + /i,/.(/i))exp{i(z(a;,+cx))-z(x,-oo))} 

= -2(1 + Hih)) exp{i(z(0, +oo) - z(0, -oo))}. 

>V(^+„^+)(x) = -2. 



Moreover, as in (15.61) 

Then, we obtain (see (|4.13l) ) 
(6.3) 



t{E, h) = {l + H{h)) exp{-(z(0, +oo) - z(0, -cx)))}, 

h 
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with 

(6.4) z{0,±oo)=i ^-^ dt. 

J ±00 



Since (j){h) is a classical analytic symbol of non-negative order and using (I2.13P we have: 

Sll 



= (1 + hm) exp{i(z(0, +00) - ziO, -00))}. 



Computation of W{lu^^^^,uj^) : As in (|6.2I 



W('^out'^in)W = 0{h)exp{--{z{x,+oo) + z{x,-oo))}. 
Using that the square root in z{x, +00) and z{x, —00) have the same determination, we have 

= exp{-i(z(0, +00) + z(0, -00))} exp{-|(z(x, 0))}, 



h 



h 



where z{0, ±00) € iM. is defined in (|6.4p and z{x, 0) = i y — ^^^^-^dt. Since the 

Wronskians are independent on x, we estimate the term z{x,0) for x = —iy, < y <^ 1. 

Here, we have z{x,0) = z{-iy,0) = -iy{i^/^^^^^^^) + 0{y^) = Cy + 0{y^) for C > 0. 
Thereafter, W{uj^^j^^,uj^) = 0{e~'^/^) for an other C > and then 

(6.5) r{E,h)=0{e-^/^). 

Consequently, using (I2.13p . we have, for a positive constant C, 



S21 



S12 = -= 



1{E, h) ^ > 
-r{E,h) 



t {E, h) 



This ends the proof of Theorem 12. 4[ 

7. Total reflection 

We suppose here that V satisfies assumption (A), the energy -E S II or IV and m > (see 
Fig. 7 or Fig. 8). 



E 



V + mc^ 



V — mc 




y+ — mc^ 




V + mc^ 



Fig. 7. E ell 



Fig. 8. E GIV 
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As in the Section we will work in two unbounded, simply-connected domains {E), 



n+{E), where Re {V{x) - E + mc^) < E, Re {V{x) 



mc 



> respectively and which 



coincide with S for |Re x\ sufficiently large. We will only work for E G II. The case E ^ IV 
can be treated similarly. Using Theorem 14.11 Proposition 13.11 and (14. 6p there are two Jost 
solutions in ^1^{E): 



exp{-z{x, 



-oo)} 



out 



exp{— — oo)} 



Hizix)) 
+H{z{x))-^ 

H{z{x)) 
-H{z{x))- 



(l + /i0(/i)) 
{l + h<i){h)) 



with 4>{h) a classical analytic symbol of non-negative order. The functions z{x, — oo), H{z{x)) 
are defined in (14.1(1 . (|5.3p and coincide with (15.41) . (|5.5p on the real axis. From Lemma WA\ 
there exist an exponentially decreasing Jost solution and an exponentially increasing one. As 
explained before Theorem 12.31 we exclude the increasing solutions which does not represent 
a physical state. We limit ourself to the one-dimensional space generated by the decreasing 
solution wj" which satisfies in Q.'^{E): 



(7.1) 



= exp{— z(x, +oo)} 
n 



-iH{z{x)) 
H{z{x))-^ 



from Theorem 14.11 Proposition 13.11 and (|4.6p . The functions z(x,-|-oo), H{z{x)) are defined 
in (|4.ip . ()3.5I) respectively and coincide, on the real axis, with 

^ [E - V(t)f - \/m2c4 -[E- y+)2 dl 



(7.2) 



z(x, -|-oo) 



(7.3) 



R{z{x)) 



{E - y+)2 X 



'mc2 ^E - V{x) 



E + v{xy 



We suppose that there is only one real turning point ti{E) and that it is simple. In that 
case the Stokes lines are as shown in the fig. Fig. 9. 





Fig. 9. Stokes lines for E GlI Fig. 10. Stokes lines for E GiV 

According to the definition of the Wronskian, we have 



a. 



out 



>V(a;+ u;, 



out,* 



>V(o. 
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Computation of yV{LO-^^,uj'^) : In order to calculate this Wronskian we need to extend one 
of the solutions u)^, uj'^ from its domain to the domain of the other solution, for example, we 
extend coj from n+{E) to n-{E). For that, we extend the square root in wj" to C\{— Im (z) > 
0, Re (z) = ti}. Thanks to the structure of the Stokes lines , we can find a path 71 from 
+00(1 — i6i) to —00(1 — (for 5i > 0) transverse to the Stokes lines along which we 
can extend w^. The extension of t £]ti{E), +oo[ — > \Jrn?c^ — [E — V(t)Y coincides with 
i^JjE^^V(fff^^rr^ on ] - cx),ii(^)[. On the other hand, on ] - oo,ti(S)[, li{z[x)) takes 
its values in e-^'^/^M+. If we denote by lo^'^ the extension of wj" along 71, we have 

u;+'^ = exp{i(-z(x,ti(i?))-z(ti(£;),+oo))}( -J/^^-^^^^ ) (1 + 
with z(ti(S),+cx)), }i{z{x)) defined in jlSl), and 



c 



(7.4) zixM^E)) = - I {{E -V{t)f -m^c^)^dt. 



On ] —oo,ti{E)[ the functions {{E — V{t))'^ —m?c^)2 and H{z{x)) are in M+. Then, we have 
W«, io+) = 2e'^l\l + H{h)) exp{i(z(ti(ii;), -cx)) - z{h{E),+^))], 

where, 2:(ti,+oo), z{ti{E),—(X)) are defined respectively in (|7.2I) and (|5.4|) . 

Computation o/yV(a;j",a;~^^).■ As in the previous paragraph we extend wj" and the square 
roots written there from Q+(£^) to Q.~{E) C C \ {Im [z) > 0, Re (z) = ti}. We can also find 
a path 72 from +oo(l — 162) to — oo(l + 182) (for 62 > 0) transverse to the Stokes lines along 

I +2 -I- 

which we can extend loJ. If we denote by oj^ ' the extension of along 72, we have 



exp{-(+.(x,ti(i?)) - .(ti(£;),+oo))} ^ 

with z{x,ti{E)), z{ti{E),+oo) and -H'(2:(x)) defined respectively in (f7^ . (f7?2]) and ([53]). On 
] — oo,ti(i?)[ the quantities {{E — V{t))'^ — m^c^)^ and H{z{x)) are in M+. Then, we have 

W{u;j,iv-^,) = -2e-'^/\l + h4>{h))e^p{^{-z{ti{E),-oo) - z{ti,+oo))}, 

where, 2:(ti,+oo), z{ti{E), —00) are defined respectively in (|7.2|) and (|5.4|) . 

Computation of W{lj^,uj~^^): Since the two solutions are defined in 17" (i?), we 

compute the Wronskian between these solutions as in (|5.6p and obtain 

W(a;r,<,J = -2. 

Then, we have 

«out = = -^(1 + ^'^(^)) eMl<h{E), -00)} 

/^J" = =e^"/'(l + /^'/'(/^))exp{l(z(ti(g),-oo)+z(ti,+oo))}, 
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with 



i rME) 

z{ti{E),-oo) = - ViE - V{tW - m^c^ " Vi^ - V-)^ - m^c^ dt 

^ J —oo 



+ -^/{E-V-Y -m?c^ti{E) 



z{ti{E), +00) = - yjm'^c^ -{E- V{t)y - y/m^c^ - {E - V+Y dt 

C J+00 

+ -^^/m?c^-{E-V+Y ti{E). 
This ends the proof of Theorem [ 



8. Zero mass case 

We suppose that m = 0, E GjV" , V'^[ and V satisfies assumption (A), (see Fig. 11) 



E 



Vz 




Fig. 11. The double turning point to 
In this section we will work in two unbounded, simply-connected domains Q,q{E), Qq{E), 
where Re V{x) < E, Re V{x) > E respectively and which coincide with S for |Re x\ suf- 
ficiently large. Repeating the constructions of the solutions (see Sections [H H]) and using 
Proposition 13. II for m = 0, there are two Jost solutions in Q^{E): 



1 



exp{—zo{x, ±00)} 



h 

-±i 



[l + h^{h)) 



~exp{^(y±-i?M 



1 



as X 



±00, 



out -exp{-^zo(a^,±oo)} 



1 

±1 



{l + hcl){h)) 



~exp{^(F± -^) x} 
he 



±1 



as X 



±00, 



with (p{h) a classical analytic symbol of non-negative order. The functions zo(x,±oo) are 
defined by 



zo{x, ±00) 



c J± 



{V{t) - V^)dt ± -{V^ - E)x. 



We suppose that there is only a simple zero of V{x) — E (see Fig. 11). Recall that the 
turning points and the Stokes lines are also defined for m = (see Definition 13.11 Definition 
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I3.2p . In our setting, there exists only a double turning point tQ{E) and the Stokes lines are 
described in Fig. 12. 




Fig. 12. The Stokes lines and the double turning point 

In this section we use the definition of the scattering and transfer matrix of Section [2] with 
the incoming and outgoing Jost solutions defined above. As in Section IH it is sufficient to 
calculate the two terms r{E,h), t{E,h) (see (|4.13ll . (|4.14l) ) to obtain the matrix S. 

Computation ofW{u!^^^^,LO^): Since the two solutions uj^ are defined in we com- 
pute this Wronskian in Qq as (|5.6p and we obtain 



out ' ^in 



-2. 



Computation ofW{uj^, io^): The two solutions Wj", io^ are defined in Qq{E), Q'^{E) respec- 
tively. Since the Wronkian W{u!^,u}^){x) is independent on x (see Remark [3.11) . we compute 
it in Qq{E) for example. For that, we extend io^ from Qq{E) to Qq{E). Using (18. ip we get 

W(Wi-,tj+) = -2{l + h4>{h))exp(^{zix,-oo) + zix,+^)) 



Then, we have 

t{E,h) = 

(8.2) 

From (I213D . it follows 



-2(1 + h4>{h)) exp ( - {z{to{E), -oo) + z{to{E), +oo)) ) . 



W«t,<) 
(1 + h(j){h))exp 

1 



(z(to(^),-oo) +z(to(£^),+oo)) 



Sll 



t {E, h) 



(l + /i0o(/i))exp{-To(i?)}, 



with 



1 / [ME] 
To{E) = -[ ( 



{V{t) - V+)dt 



ME) 



{V{t)-V-)dt + to{E){V+ -V-) 



From the form of the Stokes lines (see Fig. 12) there exist no path transverse to the Stokes 
lines along which we calculate the Wronkian >V(w~^^, wj^). More precisely, the WKB method 



24 ABDALLAH KHOCHMAN 

does not give the asymptotic behavior of the Wronkian W{uj^^^^,io^) as well as the term 

Nevertheless, (18. 2p together with the relation (|2.12|) . implies: 
(8.3) r{E,h)=0{h). 
Consequently, using (|2.13p . we have 

= WJ^) = 

-r{E,h) 



This end the proof of Theorem 12. 2[ 

Appendix A. Spectrum of the Dirac operator 

Proposition A.l. Suppose that V{x) is a application with values in the space of Hermit- 
ian 2— matrix. Moreover we assume that 

\\V{x) -V^hW ^0, as x^±oo. 

Then the operator H = Hq + V is a selfadjoint operator on D{Hq) and 

(A.l) CTessiH) = ] - oo, -mc^ + V^]U[mc^ + V~,+oo[. 

Proof. In order to prove this proposition, we first calculate the essential spectrum of Hq + W, 
where is a potential with W{x) = V^l2 for itx > i? > 0. From Lemma 5.1 of [T9j, we 
know the following inclusion: 

(A.2) aess{Ho + W) C] -oo, -mc^ + y+] U [mc^ + ,+oo[. 

Let us now prove the second inclusion. We denote 

I"*" :=] — oo, —mc? + V^"*"], and I~ := [mc? + F^, +oo[. 
For E e I^, we consider the sequence 



/± = exp{T-<I>(S-y±)x}( )x(±x/n) for n e N, 



with ^{E) = sgniE)y/E^ -m^c"^, A{E) = {J and sgn{E) = ^ for ^ [-mc^,mc^]. 

The function x ^ C'o°(I^) is such that x{x) = lif2<j;<3 and x{x) = if x < 1 and x > 4. 
The normed sequence ( n;.",, )ngN has no convergent subsequence and satisfies 

Win II 

{H^ + W -E)-^^Q, asn^+oo. 

Win II 

From the Weyl criterion, we deduce C aess{Ho + W). Consequently 

aess{Ho + W)=\- oo, -mc^ + U [mc^ + +oo[. 
Finally, using Weyl's theorem we obtain 

Cess 

and the proposition holds. □ 
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